Hartree-Fock Theory of Hole Stripe States 
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We report on Hartree-Fock theory results for stripe states of two-dimensional hole systems in 
quantum wells grown on GaAs (311)A substrates. We find that the stripe orientation energy has a 
rich dependence on hole density, and on in-plane field magnitude and orientation. Unlike the electron 
case, the orientation energy is non-zero for zero in-plane field, and the ground state orientation can 
be either parallel or perpendicular to a finite in-plane field. We predict an orientation reversal 
transition in in-plane fields applied along the [233] direction. 
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Because of the macroscopic dengeneracy of Landau lev- 
els, the physics of two-dimensional (2D) electron systems 
in strong external fields has been a fertile area for many- 
particle physics. Recently |l]-|j the emergence of strongly 
anisotropic transport properties at low temperatures has 
been interpreted as evidence for the occurrence of the 
unidirectional charge-density-wave stripe states predicted 
by Hartree-Fock theory 0. For conduction band Lan- 
dau levels with orbital kinetic energy index n > 2 (filling 
factor v > 4), the putative stripe states occur instead 
of the strongly correlated fluid states responsible for the 
quantum Hall effect ||. Although Hartree-Fock theory 
provides a clear motivation for stripe states || , it cannot 
reliably predict the nature of the ground state because 
the energetic competition with fluid states is delicate . 
Moreover, the transport properties of stripe states cannot 
be explained by Hartree-Fock theory, although they are 
consistent Sj with theories |p|dTl[| of quantum-fluctuating 
stripes. For these reasons, the ability of Hartree-Fock 
theory to predict jl2|;[l3) the low resistance (parallel to 
stripe) direction in an in-plane magnetic field has played 
an essential role in establishing the stripe-state expla- 
nation of N > 2 anisotropic transport. A recent study 
Jl4j in which a reorientation transition in a wide quan- 
tum well sample is explained by Hartree-Fock theory is 
especially convincing in this respect. 

The present work is motivated by the discovery Jl5| ] 
of anisotropic transport in 2D hole systems grown on 
GaAs (311)A substrates. In this case, anisotropic trans- 
port already occurs for v ~ 5/2, demonstrating that 
there are important differences between the electron and 
hole cases. The change is not unexpected, given the 
anisotropy of 2D hole band structure. We have general- 
ized the Hartree-Fock theory of stripe states to the case of 
valence bands described by a the single-particle Luttinger 
Hamiltonian pq| . We find that strong orbital-quantum- 
number mixing leads to anisotropic effective electron- 
electron interactions and to a dependence of stripe en- 
ergy on orientation even in the absence of an in-plane 



magnetic field. This property favors the formation of 
a stripe state, consistent with experiment. The ground 
state orientation is not in general either parallel or per- 
pendicular to the direction of a finite in-plane field when 
one is present. 

To describe 2D hole gases grown along the [311] di- 
rection, it is convenient to choose a Cartesian coordi- 
nate system with [Oil], [233] and [311] direction axes. 
The cubic coordinates in terms of which the Luttinger 
Hamiltonian is usually expressed are related to these co- 
ordinates by k a — J2i u aikii a = x,y,z; i = 1,2,3 in 
both direct and reciprocal space. Here the u a i are direc- 
tion cosines. Each element of the Luttinger Hamiltonian 
is a quadratic form in the fcj's. The two-dimensional 
hole gas is created by a GaAs (narrow gap) quantum 
well flanked by the AlGaAs (wide gap) barriers. Since 
the barrier lies in the range 100 — 400 meV and the en- 
ergy scale of interest is ~ 10 meV, we take the barrier 
to be infinite. Single particle eigenspinors of the quan- 
tum well Luttinger Hamiltonian can be expanded in the 
form tp^(x) = E lQ c ^(fc) e ik -P(i(x 3 )x a (x). Here the 2D 
wavevector k — (k±, ^2) is a good quantum number, and 
(1(^3) oc sin(?7ra;3/6) where b is the well width. Bloch 
functions Xa (%) are chosen such that when k\ and k2 are 
set to zero the Luttinger Hamiltonian takes a disgonal- 
ized form. 

For the GaAs valence bands we have used the Luttinger 
model parameters 71 = 6.85, 72 = 2.1, and 73 = 2.9 
and retained the first 20 subbands. A typical 2D band 
structure is illustrated in Fig. |l| by plotting constant en- 
ergy contours for the lowest energy subband. Since our 
confinement potential has inversion symmetry, each en- 
ergy band is doubly degenerate. For wider quantum 
wells, subband spacing is reduced and subband mixing 
is strengthened. Since cubic systems are not invariant 
under 90° rotations about the [311] direction, the 2D 
bands have lower than square symmetry. The constant 
energy contours are elongated in directions with larger 
effective mass. For the lowest subband, the Fermi sur- 
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face anisotropy is relatively weak near the zone center 
(r point), but gets stronger at larger \k\. The effective 
mass is heavier (lighter) along k 2 than along k\ direction 
close to (far away from) the zone center, and is small- 
est along the direction rotated from k\ by 45° at larger 
\k\. As we shall discuss later, the Fermi surface topology 
is manifested in anisotropic effective electron-electron in- 
teractions and ultimately in orientation-dependent stripe 
state energies. 

When a magnetic field is applied, the subband spec- 
trum consists of macroscopically degenerate Landau lev- 
els whose energies may be evaluated following famil- 
iar lines. In the Luttinger Hamiltonian, k\ and k 2 
are replaced by raising and lowering ladder operators, 
ki = i(-a + a))/\/2l and k 2 = (a + <J)/y/2t, with 
I = (hc/eB) 1 / 2 and a Zeeman term is added to the 
Luttinger Hamiltonian Hz = —k/J-bJ ■ B. Here 

fiB is the Bohr magneton and n = 1.2 is an addi- 
tional Luttinger model parameter. The envelope func- 
tion eigenspinors can be expanded in the form ipNx(x) = 



bnx(p)(i(x 3 )Xa(x) where 4>nx(p) is one of the 



parabolic band Landau gauge wavefunction generated 
by the ladder operator algebra, and the Hamiltonian 
matrix is independent of the guiding center label X. 
The eigenstates of the Hamiltonian are strong mixture 
of Landau level n and subband i indices. In diagonal- 
izing the finite-field Luttinger Hamiltonian, three sub- 
bands and 30 Landau levels were retained. This pro- 
cedure is readily generalized to allow for a magnetic 
field component perpendicular to the growth direction, 
B = B(xs + tan 9[x\ cos ip + £2 simp]). In this case, 
k\ — ► k\ + tan 9 sin px 3 /i 2 and k 2 — > k% — tan 9 cos ipx 3 / £ 2 . 

Fig. [2] illustrates our results for the magnetic field de- 
pendence of the Landau level energies for a typical quan- 
tum well width and no in-plane field. The Landau lev- 
els show very strong nonlinear dispersion with magnetic 
field. The energy levels are unevenly spaced and the ap- 
parently crossed levels are split due to the lack of parity 
under the inversion symmetry operation. Both Landau- 
level and subband mixing are stronger at larger quantum 
well widths. 

To model the stripe state seen in Rcf. at v — 5/2, 
we consider interacting electrons in the third level of Fig. 
2at5« 2.5 Tesla (n h = 1.5 x 10 11 cm" 2 ). The distorted 
semiclassical cyclotron orbit of a Fermi energy electron 
with this density is illustrated in Fig.l for the case of 
b = 250A. Semiclassical orbit distortions translate quan- 
tum mechanically into mixing of parabolic band Landau 
levels. For interactions within a Landau level this mixing 
is described exactly [Q simply by replacing the Coulomb 
interaction by 
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Note that V e s(q) is always repulsive. The parabolic-band 
plane- wave matrix elements, F nn /(q), play a crucial role 
since they are dependent on the 2D angular coordinate 
of q, 4> q : F nin MF n2 n 3 (-q) - e^+ n ^~ n ^ . It 
is precisely this effect which gives rise to orientation de- 
pendence of the stripe state energy. In an electron gas, 
rii = "-2 = "-3 = ^4 = N and the stripe state has no 
orientation dependence. 

In the Hartree-Fock approximation, the energy per 
electron for the stripe state of a half-filled Landau level 
is M 
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where 
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e is the direction perpendicular to the stripes, v* the 
filling fraction at upper Landau level, and 
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Kff(g) 
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where e(q) is the dielectric constant arising from polar- 



The evaluation of U(q) is simplified by Fourier expand- 
ing the angle dependence of V c g (q) and taking advantage 
of inversion symmetry in the 2D plane. The isotropic 
Fourier component is most dominant near q — 0. Higher 
Fourier components vanish at q = 0, and oscillate be- 
tween positive and negative values as a function of q. 
When summed, the net result is a strong suppression of 
V c g(q) at ql > 1.5, which helps to make U(q) more neg- 
ative (see Figure 3) and favors stripe states. 

The results we obtain for the orientation dependence 
of the cohesive energy in a 250A quantum well are pre- 
sented in Fig. H The ground state stripe orientation 
is tilted from the [233] axis by <f> B t ~ 23°. We studied 
the dependence of the ground state <j> on well width, b, 
finding that for b — 100A and less </> s t = 0°. On the 
other hand, for b = 150A and b = 200A Q, we find 
(f> s t w 15° and (f> s t « 22°, respectively. The difference 
in cohesive energy per electron between <f> = and s t, 
\E coh (0°) - E coh (<t> st )\, is small: for b = 150, 200, 250A it 
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is, respectively, 0.7,2.8,4.5% of the maximum difference 
|^oofc(90°) - E coh (<f> at )\. For larger values of b = 300A, 
350A and 400A, we find <p st « 45°. The ground state <f> at 
changes gradually from 0° to 45° with increasing quan- 
tum well width. Note that the dependence of the ground 
state stripe direction on well width tracks the well-width 
dependence of the Fermi surface topology. 

We believe it is band anisotropy, and the orientation 
dependence of stripe-state energy it leads to, that is 
primarily responsible for the occurrence of hole stripe 
states at v = 5/2. In 2D electron gases, stripe states 
appear [QJ^] at v = 5/2 only when anisotropy is im- 
posed by adding an in-plane magnetic field. Compar- 
ing experiment 0] and in-plane field calculations p2[ , 
we find that conduction band stripe states appear at 
v = 5/2 when the intrinsic orientation energy per elec- 
tron is ~ 0.0002e 2 /eo^, comparable to the valence band 
v = 5/2 orientation energy per hole at perpendicular 
fields illustrated in Fig. ||. 

The valence band orientation energy is also sensitive 
to in-plane field as illustrated in Fig. || for B « 2.5 Tesla. 
For quantum wells of widths b = 100 — 250 A, increasing 
in-plane field along the [011] axis (Fig. ||(a): f — 0°) 
tends to orient the easy axis of striped CDW closer to 
[233]. On the other hand, when in-plane field is applied 
along [233] axis (Fig. |(b): tp = 90°) we find that for 
a sufficiently large well widths and field tilt angles, the 
ground state orientation can change from near [233] to 
[011]. 

In comparing with experiment, it is necessary to use 
the appropriate sample width and 2D hole density and 
to account for the corrugations which occur jlTj at 
GaAs/AlGaAs(311)A interfaces. For example, we have 
investigated the hole density dependence of the orienta- 
tion energy and found that there is more anisotropy at 
higher hole densities as expected from Fig. [j]. The conse- 
quences of corrugation for the stripe states are difficult to 
quantify, partly because the stripe state period does not 
match the period of corrugation potential. Nevertheless 
it is likely that this morphological feature is partly re- 
sponsible Ea] for the factor of ~ 2 mobility anisotropy at 
zero field and that it will tend to align the stripes. Tak- 
ing account of the corrugations, our finding for the 200 A 
quantum well case ]15|] , is consistent with experimental 
finding that stripes are aligned along [233]. It is reason- 
able to conclude that the corrugation contribution to the 
orientation energy overcomes the small intrinsic value of 
\E C oh{4> = 0°) — E co h{4> s t)\- Intrinsic orientation ener- 
gies can, however, be altered by tilted fields. For fields 
along the [011] direction, E coh (<f> = 9O o )-£ co/l (0 = 0°) is 
increased and we would expect little experimental conse- 
quence. For fields along the [233] direction, on the other 
hand, we predict a dramatic reversal of stripe orienta- 
tions at tilt angle 9 ~ 60°. At this angle, the intrinsic 
band-structure effects addressed here favor [011] oriented 
stripes by - 0.0008e 2 /e ^ per electron. This number 



should be compared with the ~ 0.0001e 2 /eo^ orienta- 
tion energy produced by a 25° field tilt in n = 2 
conduction band Landau levels. In the [100] growth con- 
duction band case, stripes orient along [110] directions 
for perpendicular fields, likely because of MBE growth 
stabilities analogous to the corrugations discussed above. 
The anisotropy energy produced by a ~ 25° field tilt, in 
the n — 2 conduction band case is sufficient to overcome 
this extrinsic anisotropy and reorient the stripes. Based 
on these comparisons, we conclude that tilted field ef- 
fects can also overcome extrinsic anisotropy sources in 
the valence band case. 

In summary, valence band stripe states have a finite 
orientation energy even without an in-plane field, favor- 
ing their occurrence at v = 5/2. For quantum well widths 
b = 100 — 250A, intrinsic band effects yield a ground state 
orientation close to [233] and corrugation effects are likely 
sufficient to produce the [233] orientation seen in exper- 
iment. We predict that for typical quantum well widths, 
an in-plane field applied along the [233] direction will 
induce a stripe reorientation transition. 
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FIG. 2. Landau level dispersion with magnetic field. With 
the inclusion of the Zeeman term (k = 1.2), the spin degen- 
eracy is lifted. 
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FIG. 1. Fermi surface topology for the well width 250A. 
The constant energy contour lines (solid lines) are drawn for 
the lowest subband. fci and fa correspond to [Oil] and [233] 
directions, respectively, (for illustration, rotated images of 
solid lines by 90° are displayed as dashed lines.) The energies 
for contours are 0.61,1.1,1.71,2.31,2.92, and 3.53meV from 
inside to outside. Sample with = 1.5 x 10 11 cm -2 corre- 
sponds to the third contour (Ep = 1.71meV). 
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FIG. 3. Comparison of effective Coulomb potentials, (a) 
The Hartree potentials along [Oil] are displayed for the 
quantum well of width 250 A. Solid line corresponds to 
no-Landau-level-mixing case with screening. The upper 
(lower) dashed line represents unscreened (screened) Hartree 
potential for 2D holes at v = 5/2. (b) Hartree-Fock (HF) 
potentials, (c) Magnified view of HF potentials near the min- 
imum where the angle dependence is significant. In (b) and 
(c), solid line: without Landau-level-mixing; long dashed line: 
along [Oil]; short dashed line: along 45° from x\ axis; dotted 
line: along [233]. 
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FIG. 4. The cohesive energy at v — 5/2 is drawn as a 
function of <f> (the angle between [233] and the easy axis of 
the striped CDW) with varying in-plane magnetic field for 
the confinement potential width 250A. Solid lines correspond 
to the case without in-plane fields. 9 defines the direction of 
tilted magnetic fields from the [311] axis. The in-plane fields 
are pointing along [Oil] and [233] in (a) and (b), respectively. 
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